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Abstract 

We derive the spectra of D — 2, SU(3) supersymmetric Yang-Mills 
quantum mechanics in all fermionic sectors. Moreover, we provide exact 
expressions for the corresponding eigenvectors in the sectors with none 
and one fermionic quantum. We also generalize our results obtained in a 
cut Fock space to the infinite cut-off limit. 

1 Motivations 

Gauge theories are believed to describe all forces in Nature. One of the means 
to investigate the main characteristics of these theories is the study of their 
reduced to one point in space versions, called Yang-Mills quantum mechanics 
(YMQM) [T] [2] |3 [4 . From this viewpoint, one can calculate the zero- volume 
limit of the spectrum of the theory, then find small volume corrections if needed. 
Such approach yields the possibility to analyze the gauge dynamics which, even 
in the case of the simplest, D = 2, of these reduced systems, turns out to be 
non-trivial due to the singlet constraint. If one is interested in supersymmetry, 
it can be incorporated in the Yang-Mills quantum mechanics as well, giving rise 
to supersymmetric Yang-Mills quantum mechanics (SYMQM). 

A different motivation for studying such systems comes from the work by 
Hoppe [5] and de Witt, Hoppe and Nicolai [6], who showed that YMQM (SYM- 
QM) describe the regularized dynamics of a relativistic quantum membrane 
(supermembrane). Hence, the eigenstates of YMQM (SYMQM) turn out to be 
the wavefunctions of quantum membrane (supermembrane). However, SYMQM 
in higher dimensional spaces, where the definition of a supermembrane is consis- 
tent, are interacting theories, and difficult to solve. Thus, although the D = 2 
case is unphysical, an analytic expression for its eigenstates is of interest, since 
it may provide some indications on the form of higher dimensional solutions. 
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Several years ago a program to investigate numerically the whole class of su- 
persymmetric quantum mechanical systems in a Hamiltonian formulation by a 
cut Fock space method was proposed by Wosiek [7J . Numerous systems have al- 
ready been studied and much analytical insight was inspired by these numerical 
results [S] [5] [ID] (and the references therein). 

Although the simplest case of D = 2 supersymmetric Yang-Mills quantum 
mechanics with the SU (2) symmetry group was solved by Claudson and Halpern 
long time ago [13], very few exact solutions were found for other groups. Some 
solutions were obtained in the bosonic sector [T^] [TS] for a general SU(N) group, 
however the fermionic eigenvalues and eigenvectors seem to be unknown in the 
literature. In this work we present an analytic method inspired by the numerical 
approach of Wosiek to derive exact solutions to the D = 2 SYMQM valid in 
all fermionic sectors. We demonstrate our approach by solving the nontrivial 
system with SU(3) symmetry group. We obtain the spectra in all fermionic 
sectors and, moreover, find closed expressions for the eigenvectors in the bosonic 
sector and in the sector with one fermionic quantum. Generalizations of the 
eigenvectors to sectors with more fermions are possible. 

This paper is composed as follows. We start by introducing the class of 
supersymmetric Yang-Mills quantum mechanics and presenting the particular 
model of D = 2 with the SU(3) gauge group. This is followed by a discussion 
of the cut Fock space approach in section [3] The derivation of the solutions is 
presented in subsequent sections: in section 2] the bosonic sector is treated and 
the closed expressions for the spectrum and eigenstates are calculated; in section 
[6] the spectrum as well as the eigenstates in the sector with one fermion are 
derived. For both derivations solutions of recursion relations on the components 
of energy eigenstate in the Fock basis is needed. The theorems used to solve these 
recurrences are presented collectively in the appendices. Eventually, section [5] 
contains a summary of results. 

2 D = 2, Supersymmetric Yang-Mills Quantum 
Mechanics 

Supersymmetric Yang-Mills quantum mechanics can be obtained by a dimen- 
sional reduction of a supersymmetric, D = d + 1 dimensional Yang-Mills quan- 
tum field theory to one point in d-dimensional space. The initial local gauge 
symmetry of field theory is thus reduced to a global symmetry of the quantum 
mechanical system. The simplest ones among such systems are those obtained 
by reduction of the J\f = 1 Yang-Mills gauge theory in two dimensions [14]. 
They are described by a scalar field 4>a and a complex fermion \a, where A 
labels the generators of the gauge group. The Hamiltonian is3 

H = -n A n A + igfABC^A4>B^c, (1) 

A sum over repeated indices is assumed. 
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with the algebra of operators given by 

[4>A: ttb] = iSa.b, {Aa, A s } = 5 a ,b- (2) 
The generator of the gauge transformations is 

Ga = fABc{4>B^C ~ iA B Ac)- (3) 
The supersymmetry charges are 

Q=X A ^A, Q = A^tta, (4) 

and 

{Q, Q} = katta =2H- 2g<f> A G A . (5) 

The fermionic term of the Hamiltonian is proportional to the generator of the 
gauge symmetry and hence it is supposed to vanish on any physical (gauge 
invariant) state due to the Gauss law. Hence, the Hamiltonian is simply 

H = \rn A -KA- (6) 
Subsequently, we introduce bosonic creation and annihilation operators, 

a A = -y=((/)A + iTTA), 0L A = -j=(<f>A-WA), (7) 

and express them, as well as the fermionic creation and annihilation operators, 
in the matrix notation, 

dij = a A T^, a\j = a A Tjj, /y = f A T^, j\- = f A T^, (8) 

where are the generators of the SU (N) group in the fundamental represen- 
tation. Hence, all operators become operator valued N x TV matrices. Such a 
notation has a very practical feature, namely, any gauge invariant operator can 
be simply written as a trace of a product of appropriate operator valued matri- 
ces [16) . In situations when it is self-evident we will use a simplified notation 
for the trace of any such matrix, namely, tr(O) = (O). We get 

Qt = -L( t r(/t t)_tr(/to)), and Q = -L(tr (/at) - tr (fa)), (9) 

and 

iV 2 — 1 1 1 
H = tr (a f a) H -tr (aJaJ) - -tr (aa). (10) 

Therefore, in the case of the 577(3) group we have, 

H = tt (a f a) + 2 - itr (a f a f ) - itr (aa). (11) 

Obviously, the Hamiltonian eg. conserves the fermionic occupation number. 
Hence, we can analyze its spectrum separately in each subspace of the Hilbert 
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np = 

(a* a*) 
(gtatgt) 



Table 1: SU(3) bosonic bricks. 

space with fixed fermionic occupation number. There are 9 such fermionic sec- 
tors because of the Pauli exclusion principle. This Hamiltonian possess also 
another symmetry, namely, the particle- hole symmetry [12] 1131 , which can be 
thought of as a quantum mechanical precursor of the charge conjugation sym- 
metry in quantum field theory. This symmetry can be observed as a perfect 
matching of the energy eigenvalues in the sectors with n f and 8 — rip fermionic 
quanted. Thus, it is sufficient to consider only the sectors with np = 0, ... ,4 
fermionic quanta. 

3 The cut Fock basis method 

The cut Fock space approach turned out to be very useful in the studies of 
gauge systems for several reasons. First of all, it is a fully non-perturbative tool. 
Second, it treats bosons and fermions on an equal footing, therefore calculations 
can be extended to all fermionic sectors without difficulties. Finally, it can be 
generalized, to handle SU(N) gauge groups with N > 2, as well as systems 
defined in spaces of various dimensionality. 

The analytic approach used in this work is inspired by this numerical treat- 
ment. Hence, we now briefly summarize the numerical approach in order to 
introduce the basic notions needed in the following parts. For a more extensive 
discussion see [7] [TD] [T3] [T7] [TS] . 

The basic ingredient of the numerical approach is a systematic construction 
of the Fock basis using the concepts of elementary bosonic bricks and composite 
fermionic bricks [H] (see tables Q] and [2]) . Elementary bricks are linearly inde- 
pendent single trace operators, composed uniquely of creation operators, which 
cannot be reduced by the Cayley-Hamilton theorem. The set of states obtained 
by acting with all possible products of powers of bosonic elementary bricks on 
the Fock vacuum spans the bosonic sector of the Hilbert space [12] . As far as 
the fermionic sectors are concerned, apart of the bosonic bricks one has to use 
fermionic bricks which need not to be single traces operators |18[ . 

The spectra of the Hamiltonian are obtained by diagonalizing numerically 
the Hamiltonian matrix evaluated in the cut basis. Therefore, one introduces a 
cut-off, N cut , which limits the maximal number of quanta contained in the basis 
states. Once the cut-off is set the Hamiltonian matrix becomes finite and it is 
possible to evaluate all its elements. Subsequently, we can diagonalize numeri- 

2 This symmetry is present in D = 2. SYMQM models with SU(N) gauge group for any 

N. 
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np = 1 


rip = 2 


np = 3 


n F = 4 


(/ T « T ) 
(/t t t) 


(/W) 

(/t/tfltot) 

(/t a t a t/t a t) 

(/tot)(/t a t a t) 


(/W) 
(/ f /W) 
(/WW) 

Cf^X/W*) 
(/t a t/t/t a t a t) 

(/tat)(/Wat) 
(/t a tat)(/t/t a t) 
(/t a t a t)(/t / t a t a t) 


(/WW) 

(/^X/W 1 ) 

(/ , / t / , / t « , « t ) 

(/W)(/W f ) 
(/tot)( t/t/t/t) 

(/t a tat)(/t/t/t a t) 

(/t/t a t)(/t/t a t a t) 

(/tat)(/t a t a t)(/t/t a t) 
(/W f X/WW) 



Table 2: SU(3) fermionic bricks in the sectors with Uf = 1, . . . ,4 fermions. 



cally such matrix, as well as any other matrices corresponding to other physical 
observables. The physically reliable results are those obtained in the infinite 
cut-off limit, called also the continuum limit 17; . The procedure requires to 
perform several calculations with increasing cut-off and to extrapolate the re- 
sults to the continuum limit. It was shown that the eigenenergies corresponding 
to localized and nonlocalized states behave differently [H] [20] [S| ■ The former 
ones converge rapidly to their physical limit, whereas the latter ones slowly fall 
to zero. 

Our analytic approach uses the cut Fock basis since it provides a system- 
atic control of the Hilbert space. However, instead of evaluating the matrix 
elements of the Hamiltonian operator, one transforms the eigenvalue problem 
into a problem of finding a solution of some recursion relation. 

The analytic formulae derived in this article are particularly useful since 
they describe both the approximate numerical results and the exact continuum 
solutions. This is possible because they are parameterized by the cut-off N cut . 
In the limit of infinite N cut our formulae give the exact solutions, whereas for 
any finite cut-off they can be directly compared with numerical results. 

4 Solutions in the bosonic sector 

In this section we present the derivation of the exact bosonic solutions of the 
SU (3) SYMQM model. We use theorems presented in appendix [B] to solve 
the recursion relations obtained from the eigenvalue problem. We classify the 
solutions into sets and describe their general structure. For the simplest cases 
we demonstrate their orthogonality and completeness. 
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4.1 Recurrence relation 



A general state \E) from the bosonic sector can be expanded in the basis as 

\ E ) = a J> k (aVWaVflO). (12) 

2j+3k<N cui 

It was shown [T^] by explicit construction and by the character method, that 
the set of states {(a^a^y(a^a^a^) k \0)} containing less then N cut quanta, i.e. 
2j + 3fc < N cuU is indeed a complete and linearly independent set of states. It 
spans the Hilbert space of states with less than N cut bosonic quanta. How- 
ever, the Fock basis used in ea. (IT2l is not orthogonal, which will have some 
implications for the structure of the solutions. 

The requirement that \E) is an eigenstate of the Hamiltonian to the eigenen- 
ergy E can be translated to the condition that Qj t k must obey the following 
recursion relation (for the details of the derivation see appendix \K§ , 

aj -i, k - (2j + 3fc + 4 - 2E)a j: k + (j + l)(j + 3fc + 4)ay+i, fe 

+ l(k + l){k + 2) aj -2,k+2 = 0. (13) 

o 

The coefficients aj.k with j, k < are set to 0. Note that the first three terms 
involve the same value of the second index of aj,fc, whereas the last term mixes 
the coefficients cij_k with different values of fc. However, it is clear that since 
it only mixes fc's with the same parity, we can consider separately the sets of 
coefficients cij_k with k even and with fc odd. We will use this observation while 
classifying all solutions of eq. (fT3")l in section FOl 

The mixing term is a direct consequence of the fact that the basis states 
containing the same number of quanta are not orthogonal. Hence, solving the 
recursion relation eq. (|13|) corresponds to finding an orthogonal basis spanning 
the entire Hilbert space of states with less than N cut quanta. Comparing to 
the simpler model with SU (2) gauge group, this mixing term is a new feature. 
This follows from the fact that the SU(2) basis did not contained two linearly 
independent states with equal number of quanta. For models with SU(N > 3) 
gauge groups, there will be many such mixing terms. 



4.2 Generalized Laguerre polynomials 

Before we study in detail the solutions of the recursion eq. (fl"3|) , let us briefly 
summarize some basic properties of the generalized Laguerre polynomials, which 
will be needed in the following. 

The Laguerre polynomials C^(x) are defined as the solutions of the differ- 
ential equation 

xy" + (a + 1 - x)y' + ny = 0, (14) 
and the orthogonality relation 

/•OO 

/ (x)L%(x)x a e- x dx = 6 mn . (15) 
Jo 
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They fulfill the following three-term recursion relation, 



(n + (re) - (2n + a+ 1 - ir)££(a:) + (n + l)Cf l+1 (x) = 0. (16) 

In our problem we will encounter a recursion relation for the rescaled Laguerre 
polynomials L%(x), sometimes known in the literature as the Sonine polynomi- 
als, which are defined by, 



Then, the recursion relation eq. (| 16[) becomes 

L°_i{x) - (2n+ a + 1 - ar)L£(:c) + (n + l)(n + a + l)K +1 (x) = 0. (18) 

Equation (|18p will be called the generalized Laguerre equation, when we 
write it for a single, specific value of the index n. Contrary, when we write it for 
any n and we want to mean by it the set of all generalized Laguerre equations, 
we will call it the generalized Laguerre recursion relation. 

Moreover, we call the cut generalized Laguerre recursion relation of order 
N cut , the set of generalized Laguerre equations with N cut + 1 variables denoted 
by a (x),ai(x), . . . ,a Ncut (x). 

To simplify the notation we introduce a vector notation for the components 
aj(x) related by the Laguerre three term recursion relation eq. (fTS)) . Thus we 
will write 

S"(x) • SLj(x) = dj-i(x) — (2j + a + 1 - x)aj(x) + (j + + a + l)oj-+i(o:), 
where by S"(x) and 3Lj(x) we mean the vectors 

1, -(2j + a + l-x), (j + l)(j + a + l) ), 
a j-i( x ) 

a j+1 (x) 




4.3 Families of solutions 

It turns out that the solutions of eq. (fl"3|) can be naturally classified into separate 
sets, which we call families. Families can be labeled by a single integer k. For 
a given family, n is equal to the maximal number of cubic bricks contained in 
the basis states used to construct the solutions within this family. Additionally 
we adopt a different notation for even and odd solutions. Therefore, we define: 

• a solution belongs to the set f R if aj.k = 0, k > 2k and Oj.fc ^ 0, k < 2k. 
In words, the eigenstate can be decomposed into basis states containing 
an even number of cubic bricks and the maximal number of them is 2k, 
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• a solution belongs to the set g K if aj^ = 0,k > 2k + 1 and dj t k 7^ 0, 
k < 2k + 1. In words, the eigenstate can be decomposed into basis states 
containing an odd number of cubic bricks and the maximal number of 
them is 2k + 1. 

/o is the simplest set of solutions, for which only a^o are nonzero, i.e. is only 
composed of bilinear bricks or in other words it involves only the Kronecker 
delta invariant tensor. 

We introduce now the quantities d K (N cut ) and d' K (N cut ) which denote the 
number of states in the f K and g K family respectively at a given cut-off N cut . 
d = L-^f^J + 1 is the number of states composed exclusively of the quadratic 
elementary brick (a' a'), namely, 

|0), (aV)|0), (aV) 2 |0), (a^a^-^O). 

Hence, the energy eigenstates from the family /o will correspond to do indepen- 
dent linear combinations of such states. 

Correspondingly, there will be d ~ [ jVc " 2 t ~ 3 j + 1 states composed with 
exactly one cubic elementary brick, since there are d' such states in the basis, 
namely, 

(aVa^O), (a t a t )(a t a t a t )|0), (aV) 2 (aV a f )|0), (at a t)^-i( a t a t a t)| ). 
In general, d K and d' K are given by 



d K = 



- (N cut - 6k) 



1, d' = 



i(iV cut -6«-3)J +1. (19) 



Obviously, for a given cut-off N cut there will be in total ["^f 11 ] nontrivial 
families of solutions. 



4.4 Generic solutions 

Let us now consider generic families f K and g K . At a given, finite cut-off N cut , 
the family f K consists of d K (N cut ) solutions. Similarly, the family g K consists 
of d' K (N cut ) solutions. For each solution the eigenvector and the eigenenergy 
must be specified. We describe below, separately for the f K and g K family, 
the expressions for the eigenvectors and the quantization conditions for the 
eigenenergies. The details of the derivation are described in the appendix [51 

The form of the normalization factors is motivated by the plane-wave nor- 
malization discussed in section W. 7. II 



4.4.1 Even number of cubic bricks 

The solutions in the f K family are given by 

d K -l 

\E m ,K) even = (2E m f*e- E ™ L 6 n K+3 (2E m )(\n,2 K ) + 

K 

+ E r «-V*l n + 3 *' 2 ' s - 2 *>)' ( 2 °) 
t=i 

where the energies E m are given by the quantization condition 

Lf+\2E m ) = 0, l<m<d K . (21) 

T(x) denotes the Gamma function, T(x) = (x — 1)! for x integer and r^f" i is 
given by 

-pe^en / c\A\t—K {k ~h t)l fr%0\ 

1 24) («-t)!(2t)!' (22) 

By introducing a more convenient notation 

K 

= \n, 2k) + K-i, K -t\n + 3t, 2k - 2t>, (23) 
t=i 

we can rewrite the eigenstates in a compact form as 

d K -l 

\E m , K) even = (2E m f K e- E ™ ]T L 6 n K+3 (2E m )\^2K). (24) 

n=0 

4.4.2 Odd number of cubic bricks 

The quantization condition for this family reads 

L 6 d l +6 (2E m ) = 0, 1 < m < d' K , (25) 
and the corresponding eigenstates are given by 

d'.-l 



\E m ,K) odd = (2E m f K+ h- E ™ L 6 n K+6 (2E m )(\n,2K+l) + 

K 

+ E K-i, K -t\n + 3t, 2« - 2t + 1)) , (26) 
t=i 

where 

r^. t = (-24) f ~ % ^tut^Lr (27) 
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Again, using an simplified notation, 

|n,2« + l) = \n,2n + 1) + ^r°i d 1>K _ t |n + 3t, 2k - 2t + 1) (28) 
we can write \E m , Kj odd as 

d' K -x 

\E m ,K) odd = {2E m f K+ h- E ™ J2 L 6 n K+6 (2E m )\n,2K+l). (29) 

n=0 

In the next section we will use the above results to describe the set of all 
solutions to the eigenvalue problem at finite cut-off. 

4.5 Complete solutions at finite cut-off 

At a given cut-off N cut , the set of all solutions, {J-E)}^ , of the eigenequation 
H\E) = E\E) consists of solutions belonging to different families. As noted 
in section FOl there will be representatives of [^f^J distinct families. Hence, 
{l-E)}^ contains solutions from families f K , where 1 < n < n max = \_^N cut \, 
and from families g K >, where 1 < k' < n' max = \_\(N cut — 3) J. The dependence 
on the cut-off is hidden in the integers d K (N cut ) and d' K ,(N cut ). We now provide 
a more detailed description of {l-E)}^ . 

Since we have \n,0) — |n, 0) and |n, 1) = \n, 1) there are two families of 
solutions for which there is no mixing, namely fo and go- They consist of: 

• do solutions with E m such that L^ g (2E m ) = given by (1 < m < do) 

do-l 

\E m ,0) even = e- E ™ J2 L3 n( 2E m)\n,0), (30) 

n=0 

• d' Q solutions with E m such that L^(2E m ) — given by (1 < m < d' ) 

\E m ,0) odd = (2E m )h~ E ™ L K 2E m)\n,l), (31) 

For the remaining [^Tp-J — 2 families of solutions the mixing appears. The 
mixing becomes more and more complex as we consider families with a bigger 
maximal number of cubic bricks. We present the eigenstates from the first two 
families with a simple mixing, and then we give the expressions for eigenstates 
from the most complicated families. Hence, we have 

• di solutions with E m such that L® (2E m ) = given by (1 < m < di) 

di-l 

\E m , l) even = (2E m fe~ E ™ L 9 n (2E m )]^2), (32) 
where \n, 2) = |n,2) - i|n + 3,0). 
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• d'x solution with E m such that L^?(2E m ) = are given by (1 < m < d^) 

\E m , l) odd = (2E m )§e- E ™ J2 L n( 2E m)W^j (33) 

n=0 

where |n, 3) = \n, 3) — ^ \n + 3, 1) . 

and 

• ^K ma x solutions with E m such that L^ Kmax+3 (2_E m ) = given by (1 < 

m < d Kmax ) 

i m „-i 

I E m , l^rnax) even i^^rn) " ^ ^ ^ -^n (2£>n) 2ft ma;K ) , 

(34) 

• d' K , solutions with E m such that L 6 d ^ ma:c+6 (2E m ) = given by (1 < 

m ^ d k,J 

d' K , -1 

\E m ,K' max ) odd = (2E m ) 3 ^* +ie- E ™ ]T L n K '— +6 (2E m )\n,2 K ' max + l) 

n=0 

(35) 

The structure of solutions is shown graphically in Figure [T] Coefficients a,j t k 
are represented with dots on the j — k plane. Equivalently, each dot corresponds 
to a basis state constructed with j bilinear and fc cubic bosonic bricks. Hence, 
a cut-off with a fixed number of quanta, here N cu t — 15, can be shown as the 
oblique, straight line. The states lying below and on this line are included in the 
Fock basis, whereas those lying outside are not. The solutions of the recursion 
relation eq. p^|) are represented by dotted lines. The lowest line corresponds to 
the solutions involving only bilinear bricks i.e. the solutions from the set /o- 
The dashed line represents solutions from the set fx. The two horizontal parts 
of this line denote respectively the amplitudes dj 2 and cij . The mixing of these 
amplitudes starts at the number of quanta equal to 6, i.e. both the amplitudes 
ao,2 and 03 contain 6 quanta. 

The complete spectrum, i.e. the set of all values of the E parameter, 
{i?} , for which a nonzero eigenstate exists can be written in a compact 

form by introducing a polynomial Q^f =0 (E). The roots of 0^f =o (E) corre- 
spond to the eigenenergies of the cut Hamiltonian at a given cut-off N cutl 

{ e } n ^ = {e-- e^:°(£) = o}- (36) 



11 







I I I I I I I I 

1 2 3 4 5 6 7 

3 

Figure 1: The structure of the solutions of the recursion relation eq. (p~3|) . Each 
dot represents a coefficient a J; fc with appropriate values of the j and k indices. 
The solid lines corresponds to a sample cut-off N cu t = 15, whereas the dotted 
lines denote solutions as described in the text. 



The polynomial 0^ F (E) must be equal to the product of all quantization 
conditions relevant for a given cut-off, hence 

e£C°(£)= II ft^W ( 37 ) 

Thus, the expression eg. (1371) provides a closed formula for the spectrum of 
Hamiltonian in the bosonic sector for any finite cut-off. 

4.6 Continuum solutions 

Although, for any finite N cut the analytic spectra given by eq. (I3T1) match exactly 
the results obtained by a numerical diagonalization of the Hamiltonian matrix, 
in order to retrieve a physical information one has to perform the continuum 
limit. Numerically, this conforms to repeat several computations with increas- 
ing N cu t and extrapolate the results. On the other hand, having the analytic 
solutions we can perform the continuum limit exactly. 
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Remarkably, the solutions obtained in the preceding section for a finite cut- 
off are closely related to the continuum solutions. This can be seen as follows. 
Let us consider a set of n recursion relations with mixing such as those given by 
the recurrence eq. (1131) . Let us assume that we are interested in the solutions from 
the family f K which have the quantization condition of the form L^ffi s (2E) = 
with k < n. They are obtained by setting identically to zero all aj t k = with 
k > 2k. Their mixing coefficients are given by L^f™ t . Increasing the cut- 
off is equivalent to the inclusion of new recursion relations. However, as far 
as the family of solutions f K is concerned, the only differences occur in the 
set of possible eigenenergies, since they will be now given by the condition 
^d~(N' ) = 0> with N' cut being the new cut-off. The mixing coefficients are 
given exclusively in terms of parameters describing the recursion relations which 
does not involve the value of the cut-off. Therefore, they will be unaffected by 
its change. A similar reasoning can be applied for the solutions from the families 
9k- 

From the argument above, it follows that in order to obtain the formulae for 
the continuum solutions, one can simply extend the sum over the basis states 
to infinity. Hence, for generic eigenstates from families f K and g K we have 

oo 

\E, n) even = (2Ef*e- E L 6 n K+3 (2E)\^2K) (38) 

oo 

\E, n) odd = (2Ef*+h- E ]T Ll K+6 (2E)\n,2K + l) (39) 

ra=0 

Some of the properties of these solutions are summarized in the next section. 



4.7 Properties of \E,n) even solutions 

In this section we verify the orthogonality and completeness of the solutions 
derived in the preceding section. Since we want to demonstrate the mechanisms 
that are responsible for these properties we will only deal with even solutions. 
Therefore we neglect their subscript even- The discussion for the odd solutions 
follows the same lines. 



4.7.1 Orthogonality 

The scalar product of two eigenstates of the Hamiltonian reads 
(E,a\E',/3) = e- E - E ' {2Ef a {2E'f p x 

^ ^ F(n + 6a + 4)F(m + 6/3 + 4 V ' 1 ' Hl V ' 

n=0 m— 

The general expression for the scalar products of orthogonalized basis states 
(n, 2a\m, 2/3), is not known. However, one can calculate them for few simplest 
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cases (namely, a = 0, 1, 2). We find that 



(n, 2a\m, 2/3) = 5 a p5 nm N a (n), (41) 

and 

Af°{n) 
Af 2 (n) 



r(n + l)r(n- 


h4) 


6 


* 


r(n + l)r(n- 


-10) 


3456 




r(n + l)r(n- 


-16) 



1990656 

which can be summarized in a compact form as 

Ua{n) = \^ T{n + 1)r(n + 4 + 6a) ' 
We can immediately conclude that solutions belonging to different families are 
orthogonal, even if they have the same energy. 

As far as the orthogonality of solutions within each family is concerned, we 
have 

(E,a\E,(3) = 8 af) e (4EE ) -_ ^ r(w + 6a + 4) ^ ! ( 42 ) 

m— y 

The finite sum in eg. (l4"2l can be calculated [55] and yields 

-E-E'/twnSa 1 1 K) ! 1 



<f?,a|^,/3) = 8 a0 e- E - E {AEE 



6 24 2 « T(d Q + 6a + 4) 2{E - £') 

X (£^+3(2^)£6a+3( 2S ')_£6a+3^ £ 6a+3( 2 ^/)Y (43) 

Hence, for and £" fulfilling the quantization conditions, £® a+3 (2E') — and 
C 6 d " +3 (2E) = 0, the scalar products reads 

(E,a\E>,f3) = J Q ^ EE ,l-l-4a ± 3 (2i;)£ 6a+4 (2i;) (44) 

The fact that the eigenstates are orthogonal among each other is a simple con- 
firmation that we managed to construct a set of orthogonal states in the Hilbert 
space with even number of quanta. 

The question of orthogonality of the continuum solutions is more subtle. The 
continuum limit of even solutions is given by eq. (|38p . Following the same steps 
as for calculations with finite cut-off and using the formula [55] for an infinite 
series of product of two Laguerre polynomials we obtain 



1 1 



{E,a\E',P) = S a ^-^e- E ~ E (4EE') 3a 



' 6 24 2q 



x 



x lim 



(1-z) L cxp(-z-^ —^-){AEE'z) 2 7 6q+3 ^ _ 



(45) 
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Introducing e = |(1 — z) and exploiting the asymptotic behavior of the Bcsscl 
function for large argument, namely 



I a (x) 



'2irx 



for large x, 



we get 



lim 



1 



(E,a\E',/3) =5 a 8-UEE')-^-—. s — 



e * 



Eq. (|4"Tl) reduces thanks to the well-known Dirac-delta representation 

1 



lim 



S(x) 



to the desired result 



(E,a\E',p) = S a/3 ^^(2E)- 3 S(V2E~ y/2E>) 



(46) 



(47) 



(48) 



(49) 



12 24 2a 

The norm of these eigenstates is divergent, which is in agreement with our 
expectation that our solutions in the continuum limit should behave like plane- 
waves. 



4.7.2 Completeness 



Next we show that the set of states derived above is a correct basis of the Hilbert 
space. Such set of solutions must be complete, i.e. any other solution can be 
written as a linear combination of \E, k). This can be done by showing that 
any Fock state can be expressed in terms of the energy solutions. Hence, the 
matrix representation of the transformation of the set of Fock states into the 
set of solutions \E, k), denoted by T, must be nonsingular. 

T can be written as a product of two matrices, T = T2T1, where T\ represents 
the matrix representation of the transformation { \j, k) } — > { \j, k) } , and T2 rep- 
resents the matrix representation of the transformation {|j, k)} — > {\E m ,n)}. 

It can be seen from the definition of the orthogonalized Fock states eq . (f23]) 
that the matrix T% is a triangular matrix with all diagonal elements equal to 1. 
Hence, its determinant is simply 1. 

The matrix T2 is block diagonal with each block corresponding to a given 
family of energy eigenstates \E m , k). Hence, the determinant of T2 will be given 
by a product of determinants of the transformation matrices written for each 
family independently. To calculate the latter at a finite cut-off, let us consider 
d K solutions from the f K family. By definition, they can be obtained from the 
orthogonalized basis states through the following matrix 



\Ex,a) 
\E 2 ,a) 



\ ( 



\ \E dK ,a) J 



6«+3 


6«+3 




(2£i) 
(2E 2 ) 



Lf+liZEi) \ 



V L 6 K+3 (2E d 



( |0,2k) \ 
II, 2k) 



LZ + _\{2E dK ) J \ |d K -l,2K) J 
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One can show that the corresponding determinant is the Vandermonde deter- 
minant depending on energies E m . To see this let us write the Laguerre poly- 
nomials Cl K+3 (E) as 



-6k+3 



(2E) 



1 



i!r(n + 6K + 4) 



(225)' 



n-l 

E 

fc=0 



c n , k L 6 k K+3 (2E) 



where c n ^ are some constants. Then, adding to each column an appropriate 
linear combination of preceding columns, starting by the last one, and ending 
on the second one, we get the determinant 

(d„-l)!r(6K+3+m) (2-^1 ) d " 
(d K -l)!r(6K+3+m) ( 2jB 2) 



2E X 
2E 2 

r(6 K +4) r(6 K +5) 2 ^ d « 



T(6k- 


-4) 


T(6k+5) 
1 


T(6k- 


hi) 


T(6k+5) 


1 




1 r 



• [3« 

d K -l 

(n 



fc=0 



-l)!r(6/c+3+m) 
1 

fc! F(6k + 4 + fe) 



IJ (2Ej - 2Ek) 

l<i<j<d K . 



Recall, that the energies E m are given by the quantization condition £^ K+3 {2E) = 
0. The roots of the Laguerre polynomial are simple, i.e. the set of roots does 
not contain two equal numbers. Hence, the above Vandermonde determinant is 
nonzero. 

In order to extend this result to the continuum limit it is necessary to prove 
that the Laguerre polynomials in the limit of infinite order have only simple 
roots. This can be done using the following formula )22j . 



lim n 

n— >oo 



-6k— 3 /*6/c-f-3 



6k+3 



(50) 



The Bessel function of the first kind, J$ K +3(2^), has countably many simple 
zeros. It follows that the zeros of the Laguerre polynomials in the limit of infinite 
order are also simple. Therefore, the infinite Vandermonde determinant will be 
also nonzero. 

The determinant of the T matrix is the product of determinants of T\ and T% 
matrices. Hence, we showed that it is nonzero for both finite and infinite cut-off. 
This proves that the transformation matrix T is nonsingular and therefore the 
set of solutions {|-E m , k)\ is a correct, orthogonal basis of the Hilbert space for 
the system with the SU (3) gauge symmetry. 



5 General recurrence relation in fermionic sec- 
tors 

In this section we derive the recursion relation for the components of energy 
eigenstates in the Fock basis in the fermionic sectors. We find a general recur- 
rence in a sector with tif fermions, and then specialize to the case of np = 1 in 
the following section. 
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In order to proceed we must introduce some notation for the fermionic bricks 
such as those presented in table [2j Let us denote a generic fermionic brick in 
the sector with n F fermionic quanta by C^(n B , Up, a). The index a labels all 
fermionic bricks in this sector and takes values in 1 < a < d nF , where d nF is 
the total number of fermionic bricks the n F sector, n B designs the number of 
bosonic creation operators contained in the a-th fermionic brick. 

The basis in the sector with n F fermions is composed of the following states 

\j,k,a) = C*(n%,n F ,a) (aV^VaVflO), (51) 
Hence, a general state can be decomposed as 

\ E )=J2 E *j,k\jM (52) 

a=l 2j+3k<N cut -n<^ 

Let us rewrite the eigenequation for H in terms of the coefficients a" k , namely 
(H-E)\E) = 

E E a lk ([H,tf(n%,n F ,a)] + C\n B ,n F ,a) H)\ 3 ,k) = 0. 

a=l 2j+3k<N cut -n% 

Evaluation of the commutator yields 

[H,Cft(n%,n F ,a)] =^n B tf(n B ,n F ,a) -±[(aa),C*(n B ,n F ,a)] (53) 

We must now evaluate the commutator [(aa), (n B , n Fl a)] . For a particular 
(n B , n F ,a) it will be given by a sum of n B terms, each equal to (n B , n F ,a) 
with one of the bosonic creation operators replaced by a bosonic annihilation 
operator. Let us denote these terms by G l al where the index t runs from 1 to 

[(aa),CHn B ,n F) a)] = ( 54 ) 

In order to get rid of the annihilation operators contained in G l a , we would like 
to push them over the creation operators (a) a*) 3 (a) a) a)) k , so that they hit the 
Fock vacuum. We have 

V* [G^aV^aWf] = (aV) 5 ' [G* a , (aW)*] + [G* Q , (aW) J '] (a t a t a t ) fe 

(55) 

Since, G l a contains exactly one annihilation operator, one can simplify the above 
commutators, as 

V t [Gl (oW)'-] = j (aV^ 1 [G l a , ( t t)] , (56) 
V t [&„, (aW) fc j = k ( a t a t a t)fc-i [ G * qj (oW)] . (57) 
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Therefore, 

V t [G^aV^VaVf] = k [G^,, (a^aV)] (aV^VaV)*" 1 

+ 3 [G^C^a 1 )] (aV^'-^aW)* (58) 

Eventually, it is easy to evaluate [G^, (aW)] since this simply yields back the 
G' (rig, np, a) composite fermionic brick. Thus, we obtain 

£ J2 {c\n%n F ,a)(s^+^-*l k 

a = l 2j+3k<N cut -n% 

+ l(k + l)(k + 2)aj_ 2!k+2 )\j,k) 

+ H<k(k [Gt,(at t a t)] |i,fc-l) + (aVy(atatat) fc G^|0))}=0 
t=i 

(59) 

Remarkably, we recover a Laguerre recursion relation for the a" fc coefficients 
with an index shifted by n% compared to the bosonic case. The mixing terms 
produced by [G^, (a^aW)] mix coefficients a" fc with different a. Note that 
[G^, (aWa'f)] cannot be proportional to G^ (rig,np,a) since it contains one 
additional bosonic creation operator, whereas the part of G* a which does not 
annihilate the vacuum misses one creation operator compared to C^(rig, rip, a). 
Therefore, the two last terms of the recursion relation are true mixing terms. 
We will investigate in details these recurrence relations in a particular case of 
the sector rip = 1 in the next section. 

6 Solutions in the rip = 1 sector 

In the sector np — 1 we have two fermionic bricks, namely, 

(/t a t) = Ct (1,1,1), (/W) = Gt(2,l,2). 

Hence we have, 

[M,(./V)] =(pa) = G\, 
[(aa), (/W)] = (/t t ) + (/W) = G\ + G\. 

Then, 

[G^aV^aWf] = ^(/t a V)(aV)>W) fe ~ ] 
+ j (/ t a t )(a t a t )^ 1 (a t a t a t ) fe , 

[G\,{a)(Jy{a)(Ja)) k ] = ^(/t a t)( a t a t) J +i( a t a t a t ) fe-i 
+ 3 (pa^a^(a^a^y- 1 (^a^a^ k 
= [G2,(aV)->( a t a t a t)fc] 
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(60) 
(61) 

(62) 



This gives the set of recurrences 



S f +4 ' 4* + + !)( fc + 2 ) «)-2, fc+2 + ^ a ?-i, fe +i ^ 0. 



Sf + 5 • a 2 , fe + 3 -(k + 1)(* + 2) a?_ 2>fc+2 + ^*±*> = 0- 



(63) 



The next step is to decompose these recursion relation into even and odd 
parts. 

The odd part contains equations for k with k even and a| k with k odd. 
Therefore, we set k = 2m and fc = 2m + 1 respectively in the first and second 
recursion relation of eas. (|63[) and obtain 

Sf l+4 • 4 2m + g (2m + l)(2m + 2) a)_ 2<2m+2 + o?_ 1)2m+1 = 0, 

S - m+8 ' a i,2m+i + g(2m + 2)(2m + 3) a)_ 2 ^ m+z + 3(m + 1) a\ 2m+2 = 0. 

Since, the coefficients a] 2m appears always with an even second index, and the 
coefficients Oj 2m +i appear always with an odd second index, the superscript 
becomes redundant and can be omitted. Hence, we write a] „ , — > a] m and 
4 2m+i -> a f 1+1 and get 



g 6 m+ 4 . , m + 2m+l_ + 3 + + 2m + 2 = 0j 

2 o (64) 

S j m+8 ' a f l+1 + 3 ( m + X ) fl f n+2 + |( 2m + 2 )( 2m + 3 ) a |-2 +3 = °- 

We proceed similarly for the remaining set of equations. After setting 
k = 2m + 1 and k = 2m in the first and second recursion relation of eas.(|63| 
and changing the notation aj 2k — > a^ fe and a] 2fe+1 — > a j fc+1 ) we have 

s 6m+5 . a 2m + 3(2m + 1) &2m+1 + 3 + + „ +2 = ^ 

2 § (65) 

S 6m+7 . a 2m+l + ^ + ^ ffl 2m+2 + | (2m + 2)(2m + 3) fl 2m+3 = q 

Before writing the expressions for the solutions of the above recurrence re- 
lations, we discuss the general structure of the solutions. 



6.1 Families of solutions 

In a similar manner as we did in the bosonic case, it turns out to be possible to 
classify the solutions in the np = 1 sector into families. We have just argued 
that parity is a good quantum number. Furthermore, the set of solutions of the 
same parity can be divided into two distinct sets of families. 
We adopt the following notation: 
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• a solution belongs to the family g\ if it is composed of basis states with 
an odd number of quanta and the basis state with the maximal number 
of cubic bricks equal to 2k, is proportional to the (/W) brick, 

• a solution belongs to the family g 2 if it is composed of basis states with 
an odd number of quanta and the basis state with the maximal number 
of cubic bricks equal to 2k + 1, is proportional to the (pa 1 a 1 ) brick. 

The recursion relation for the families g\ and g 2 are given by eas. (|64l) . 

• a solution belongs to the family f\ if it is composed of basis states with 
an even number of quanta and the basis state with the maximal number 
of cubic bricks equal to 2k + 1, is proportional to the (pa') brick, 

• a solution belongs to the family f 2 if it is composed of basis states with 
an even number of quanta and the basis state with the maximal number 
of cubic bricks equal to 2k, is proportional to the (pa' a') brick. 

The recursion relation for the families /* and f 2 are given by eas. (|6"5j) . 



One can easily calculate the number of solutions within each family, by 
counting the number of basis states fulfilling appropriate conditions. Hence, in 
analogy to the bosonic case we define d\{N cut ) and d 2 (N cut ) as the numbers 
of solutions in the families j\ and f 2 at cut-off N cut . Similarly, d'^(N cut ) and 
d' 2 (N cut ) denote the numbers of solutions in the families g\ and g 2 at cut-off 
N cut . They are given by 



d'Z(N cut ) 



dl(N cut ) = 



d'Z(N cut ) 



^(N cut -6K-1))\ +1, 
i(JV ctlt -6«-4))J +1, d\(N cut ) 



\{N cut -6n-b)) 



+ 1, 
(66) 

i(JV cut -6K-2))J +1. 

(67) 



Obviously, for a given cut-off N cut there will be in total |_ 
families of solutions. 



N cut 



6.2 Generic solution and the complete set of solutions at 
finite cut-off 

In this section we use theorem [2] from appendix iBl to write down the solutions 
to the recursion relations eqs. fM)) and eqs. (p5|) . We treat a generic cases of 
solutions belonging to the families g\ and g 2 . and then to the families j\ and 
f 2 . Eventually, we describe the complete set of solutions at finite cut-off. 

6.2.1 Solutions with odd number of quanta 

The quantization condition for the family g\ reads 

L% + \2E m ) = 0, l<m<4, (68) 
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and the corresponding eigenstates are given by 

d'}-l 



\E m:K ,l) odd = (2E m ) 3K+ h- E - J2 L 6 n K+i (2E m ){(pa^)\n,2K) + 

K 

+ E ( r 2S K -2p + i(/ t « t a t )h + 3p - 2, 2« - 2p + 1) + 

+ m d 2 K -2 P (P^)\n + 3p, 2« - 2p)) }. (69) 



p=i 



Similarly, the quantization condition for the family g 2 reads 

L^ +8 (2£ m ) =0, 1 < m < < 2 , (70) 
with the eigenstates if the form, 

d' 2 -i 



\E mi K,2) odd = (2E m ) 3K+ h- E ™ J2 L 6 n K+s (2E m ){(pa^a^)\n,2K+l) + 

n=0 

+ E ( r 2Si, 2K - 2p+ 2(/ t fl t )k + 3p - 1, 2k - 2p + 2) + 

+ r^ +1 , 2ti _ 2p+1 (/ t at a t)|n + 3p, 2k - 2p + 1)) }. (71) 



P =i 



rjfj are given by recursive relations in theorem [2] presented in appendix [B] 

6.2.2 Solutions with even number of quanta 

In this case, the quantization condition for the family /* reads 

L^+ 7 {2E m ) = 0, 1 < m < 4, (72) 
and the corresponding eigenstate is of the form 

di-i 



\E m ,K,l) even = (2E m ) 3K+2 e- E ™ E Ll K+7 {2E m )[{fa))\n,2K+ 1) + 

n=0 

K 

+ E ( r 2:Ti, 2K -2 P (/ t « t a t )l« + 3p + 1, 2« - 2p) + 

+ r^ e +i, 2K -2 P+ i(/ t « t )h + 3p, 2 K - 2p + 1)) }. (73) 



p=i 



The quantization condition for the family is given by 



L% +5 (2E m ) = 0, 1 < m < dl (74) 



21 



dt-1 



and the eigenstate reads 

\E m ,K,2) even = (2E m )^ +1 e- E ™ J2 L 6 n K+5 (2E m ){(pa^a^\n,2K) + 

K 

+ E ( r 2^-2p + i(/ t « t )l^ + 3p - 1, 2« - 2p + 1> + 

+ r^_ 2p (/W)|n + 3p, 2« - 2pj) }. (75) 



p=i 



Again, are as given in theorem [5J 

6.2.3 Complete set of solutions at finite cut-off 

In a similar manner as in the case of the bosonic sector the complete set of 
solutions at a given cut-off iVcutj {1^)}^ IS given by the union of all solutions 
from nonempty families. We write explicitly only the d'Q(N cut ) odd and d^{N cut ) 
even simplest eigenstates with no mixing, namely, 

do 1 -! 

\E m ,0,l) odd = (2E m )h- E ™ J2 Li(2E m )(pa^)\n,0), (76) 



n=0 



with 
and 

with 



L*n(2E m ) = 0, l<m<d' \ 

4-1 

\E m ,0,2) even = (2E m )e~ E ™ £ L 5 n (2E m )(f a^)\n, 0), 



(77) 



n=0 



L 5 d2 (2E m )=0, l<m<d 2 . 



The spectrum, i.e. the set of all values of the E parameter, {i?} , for 
which a nonzero eigenstate exists, can be written in a compact form using the 
polynomial e^f^E). The latter can be deduced to be of the form 



li(JW-i)J 



n £ S 

J=0 



LitiVou t -3i-i)J+i( 2jEJ ) 



L|(iVe« t -2)J 



n i ?i(L t -3 fc - 2) j +1 ( 2 ^).- ( 7g ) 



fc=0 



Hence, again, the expression ea.(|78|l provides a closed formula for the spectrum 
of Hamiltonian in the fermionic sector for any finite cut-off. 



22 



6.3 Continuum solutions 



The argument described in section 14.61 applies as well to the case of recursion 
relations eqs.(|63p. Therefore, one can obtain exact, continuum solutions from 
the solutions described in the preceding section, by simply extending the sums 
to infinity, i.e. solutions eqs. (|7o) and (fTTf become 

oo 

\E, 0, l) odd = (2E)h- E £ 4(2£)(/W)|n, 0), 

(79) 

\E,0,2) even = {2E)e- E Y J Ll{2E){pa^a))\n 1 0). 

71 = 

The energy quantization conditions in the continuum limit allow energies that 
lie on the whole real, positive axis. 



7 Spectra in all fermionic sectors 

In section[S]we derived a general recurrence relation valid in all fermionic sectors. 
Remarkably, it has the form of a set of generalized Laguerre recurrence relations 
coupled by a number of mixing terms. However, for any finite cut-off iV cut there 
always exists one homogeneous Laguerre equation for which the mixing vanishes 
since it involves terms with a number of quanta beyond N cut . By solving these 
homogeneous equations and then proceeding with the remaining equations (in 
analogy to the proofs of theorems [T] and [5]), it is possible to obtain solutions of all 
recursions in all fermionic sectors. In particular, the corollary[T] which justifies 
the existence of a compact polynomial (E) which roots corresponds to the 

eigenenergies of the Hamiltonian, can be extended to all sectors. 

The polynomial (E) is equal to the product of quantization conditions 

of all nonempty families of solutions at any finite cut-off N cut in a given fermionic 
sector. Each quantization condition is characterized by two integers, n and 7, 
and has the form L^{(2E) — 0. The index 7 is given by the general recursion 
relation, i.e. 7 = 3fc + 3 + n%{np) : where for each i, the index k enumerates 
the families of a specific type and is bounded by < k < [^(N cut — rig(n F ))\. 
The index n corresponds to the number of solutions within each family and it 

is equal to n = ^(iV CMt — 3fc — n^(rii?))J + I. Hence, the general polynomial 
can be written as 



\i(N cut -3k-n%(n F ))i+l 



^i t ^)=n ' n ^:: s 3rU(„ F )), +1 (^))- m 

a=l \ k=0 

Therefore, in order to be able to determine the spectra the numbers ng(np) 
must be known. In the case of the SU(3) model they can be simply read 
from the table [2] The above prescription for the spectra was crosschecked with 
independent numerical calculations up to cut-off N cut < 40, in all fermionic 
sectors < np < 4. Both computations yielded exactly the same numbers. 
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8 Conclusions 



In this article we have studied the model of D = 2, Supersymmetric Yang- 
Mills Quantum Mechanics with a SU(3) gauge group. We have proposed a 
method for solving this system. It applies as well to the bosonic sector as to 
the fermionic ones. We have derived explicit formulas for the eigenstates in the 
bosonic sector and in the sector with one fermion by exactly solving a recursion 
relation involving mixing terms. Such mixing terms are a novel, characteristic 
feature of the SU(N > 2), SYMQM models. Moreover, we obtained expressions 
giving the energy spectrum in all fermionic sectors which are parameterized 
by the cut-off. Hence, on one hand, one has an analytic confirmation of the 
correctness of numerical calculations, which are always done at finite cut-off. 
On the other hand, one can easily and in a precise way extrapolate these results 
to the continuum limit and therefore obtain expressions for the exact spectrum 
and corresponding eigenstates. 

These results can prove to be valuable in studies of several further problems. 
A generalization to solutions of the models with SU(N) gauge groups with 
N > 3 is possible not only in the bosonic sector but in fermionic sectors as well. 
The results for general N may be very useful in the investigation of the large- N 
limit. Furthermore, in the search of solutions of more physically interesting, i.e. 
D = 4 or D — 10 dimensional models, one can try to derive and solve a similar 
recursion relations. 
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A Derivation of recursion relation eq. ( |T3| ) 

We have, 




(81) 



and 



oo oo 



\E) = ]T a id \i,j) = MaVnaWnO). 



(82) 



i,j—0 i)i=0 
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One can calculate that, 

(ataXaVyVaVylO) = (i + |?) (aV) VaWlO), (83) 

(aVXaV^aVa^'lO) = (aVf+^aVa^'lQ), (84) 
(aa)(a t a t )' i (a t a t a t ) J |0) = 3.? + 3)(a t a t ) l - 1 (a t a t a t ) J |0) 

+ - l)(at a t)^ (a t a t a t )J -2| 0)j (85) 

where for the last equality we have used the following commutators, 

[(aa), (aW)*] = t(i + 3)(a t a t )^ 1 + 2i(aV) 1 - 1 (a+a), (86) 

[(aa), (aV a t)*l = |t(*-l) (aV a ty-2( a t a t)2 + 3i ( a t a t a t)*-i( a t a t a ) 5 ( 87 ) 
8 

and 

[(aVaM^aV)*] = ^(aVa^^aV) 2 (88) 

Thus, 

ff|i,J> = (i + |i + 2)|i,j)-i|< + l,i>+ 

- + 3j + 3) |* — - - l)|i + 2, j - 1) (89) 

Therefore, grouping the coefficients in front of each basis state yields the recur- 
sion relation 

ai-i d - (2i + 3j + 4 - 2£)a iJ + + + 3j + A)a l+1J + 

+ |0' + l)0' + 2)ai_2j+2 =0. (90) 

B Solutions of the cut recursion relations 

In this section we solve cut recursion relations for the generalized Laguerre 
polynomials. We construct the theorems in the order of increasing complexity. 
We start by providing the solution to a single cut recursion relation (lemma 
[I]). Next, we consider respectively, the problem of a single inhomogeneous cut 
recursion relation (lemma [5]), the problem of two coupled recursion relations 
(lemma [3]) , and the problem of a set of coupled recursion relations (theorem [T]) . 
Eventually, we solve a problem of a set of twofold coupled recursion relations 
(theorem and formulate a corollary (corollary [1]), which can be used to get 
the spectra of the system in all fermionic sectors. 
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B.l Some lemmas 

Let us start by considering the following lemma, 

Lemma 1. Consider a set of k + 1 generalized Laguerre equations for the coef- 
ficients ao, . . . , a K , 

aj-i(x) - (2j + a + 1 - x)a j {x) + (j + + a + l)a j+1 (x) = 0. (91) 

Then, there exist k + 1 nontrivial solutions, denoted by (a.,) , . . . , (oj) K , of the 
form 

(a 3 y =a Q T(a + l)Lf(xi), (92) 
where Xi are the solutions of the equation 

L a K+1 ( Xi ) = 0, i = 0,...,/t. (93) 

Proof. Let us write this cut set of generalized Laguerre equations in an explicit 
form as 

— (a + 1 — x)a + (a + l)ai = 0, 
a - (a + 3 - x)ai + 2(a + 2)a 2 = 0, 
ai — (a + 5 — x)a 2 + 3(a + 3)a 3 = 0, 
a 2 - (a + 7 - x)as + 4(a + 4)a 4 = 0, 

a K _ 4 — (a + 2k - 5 - a;)a K _ 3 + (k — 2){a + k — 2)a K _ 2 = 0, 
a K _ 3 - (a + 2k - 3 - a;)a K _ 2 + (« - l)(a + k — l)a re _i = 0, 
a K -2 - (a + 2k - 1 - a;)a K _i + k(q: + /t)a re = 0, 
a K _! — (a + 2k + 1 — x)a K = 0. 

Now, one can express a\ in terms of ao using the first equation, then express a 2 
in terms of ao using the second equation. Proceeding in this way for all but the 
last equation, enables us to rewrite the above set of equations as 

Oi 

a 2 
as 

a K -i 
a K 
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a + 1 



-ao 



a + 1 

((a + 3-a02±i=£-l) 



2(a + 2) 
a r(a + l)Lf(a;), 
a T(a + l)L2(x), 



a r(a + l)L?(x), 

a = a r(a + 1)L 2 (x), 



:a r(a + l)L«_ 2 (ar), 
:a r^ + l^^x), 
:a r(a + l)L°(x). 



The last equation reads, 

a £«_i(z) - a (a + 2k + 1 - x)L%(x) = 0. (94) 

Using eq. ([T5)) it can be transformed into 

a L« +1 (x)=0 (95) 

Thus, if x is tuned to be a root of then ao can be arbitrary and the system 

admits a nonzero solution. On the contrary, if x is not a root of L", lt then ao 
must vanish. Hence all coefficients dj vanish too. 

Eventually, the equation L" +1 (x) — has k + 1 zeros, so the system has 
k + 1 independent solutions. We denote them by 

(a j )°,...,(a j r (96) 

with 

(a 3 y =a Q T(a + l)L«(x l ). (97) 

□ 

We now prove two lemmas giving the solutions of an inhomogeneous Laguerre 
recursion relations. 

Lemma 2. Consider a set of n+1, inhomogeneous, generalized Laguerre equa- 
tions of the following form, 

aj-i(x) - (2j + a + 1 - x) aj (x) + (j + l)(j + a + l)o i+ i(x) = X^ +q (x), (98) 

with x ^ 0, a and (3 such that Mf3 — a) is a positive integer and x not being a 
root of the polynomial L" +1 (x). 

Then, there exist specific values of a, f3 and q, for which the system admits 
solutions with aj nonvanishing. 

Proof. Equations (|98j) form a set of k + 1 inhomogeneous equations with n + 1 
variables a,j(x). Its determinant is equal to L" +1 (x) and, by assumption, does 
not vanish. Therefore, the system admits a unique solution. We will construct 
it in the following way. We assume that a.j (x) is proportional to Lj_ p with some 
proportionality factor 7 ^ and p some integer to be fixed latter, 

aj (x)= 7 L?_ p (x). (99) 

Then, the recursion relation for aj (x) takes the form 

if_ p _x(x) - (2j + a + 1 - x)L^_ p (x)+ 

+ (J + l)(j + a + l)L?_ p+1 (x) + ^ +q (x) = 0. (100) 

The general recursion relation for Lj is 

LP_ x {x) -(2j + (3 + l- x)L^(x) + (j + 1)lP +1 (x) = 0. (101) 
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By assumption we have that g(/3 — a) = k with k integer, k > 0. Thus, for 
j > k we can shift the index j — >• j + k in eq. (|100|) . After setting p = k, we 
obtain, 

L$_ 1 (x)-(2j+l3 + l-x)L$(x)+ 

+ (j + k + l)(j + k + a + l)L? +1 (z) + ~L? +k+q (x) = 0, (102) 

which can be considerably reduced using eq. (|101[) . Thus, eventually, 

-i(/3 2 - a 2 )L? +1 (x) + ^ +k+q (x) = 0. (103) 

Eq. (1103|) must be satisfied for any value of the x parameter. This can only 
happen when 

q=l-k = l-±(J3-a), (104) 

7=732^- ( 105 ) 
— or 

It follows from eq. ([99]) that the coefficients a,j (x) for j < k must vanish. 

Summarizing, we showed that there exist a nonvanishing solution to the 
initial set of equations, given by 

aj (xi)=0, j = Q,...,k-l (106) 
aj(xi) = yL?_ k (xi), j = k,...,K, (107) 

where k = a). Note that this solution has no more the freedom of arbitrary 
constant factor. □ 

Lemma 3. Consider two cut, generalized Laguerre recursion relations, one of 
them supplemented with a mixing term, \i °f th £ following form, 

a i-xi x ) - (2j + a + l- x)a,j(x) + (j + + a + l)a j+1 (x) - xbj+ g (x) = 0, 

b^x(x) -(2j+p + l- x)b 3 (x) + (j + + 13 + l)b j+l {x) = 0, 

(108) 

and let us assume that — a) = k with k a positive integer and q < 0. 

Then, for some specific values of the parameter x the system admits solutions 
with both aj (x) and bj (x) nonvanishing. 

Proof. The set of equations for bj (x) can be brought the form of lemma [TJ 
Therefore, for the cut-off k, it admits k + 1 nonzero solutions. Consistency 
requires to assume that there are k + 2 — q coefficients a,j (x) . Considering the 
possible values of the x parameter we have three cases: 

• ^K+2-q( x ) 7^ an d L^ +1 (x) 7^ 0, then neither, aj(x) nor bj(x) admit a 
nonvanishing solutions. 
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• L" +2 _ q {x) — and u K+1 [x) ^ 0, then bj(x) vanish, yielding the homoge- 
neous set of equations for a,j(x). Because L" +2 —q( x ) = 0, a-j(x) will admit 
k + 1 — q nontrivial solutions of the form, 

(a 3 Y = a r(a + l)L"(x,), j = 0, . . . , k - q + 1 (109) 

with being the solution of L" +2 _ q (x) = 0. 

• ^k+2- 9 ( x ) 7^ and = 0, then bj(x) admit k + 1 nontrivial so- 
lutions. However, due to the mixing term, a,j{x) cannot vanish. As was 
shown by the preceding lemma[2j for 2 {P~ a) = 1 — q there exists a unique 
solution to the initial set of equations with dj(x) and bj(x) nonvanishing, 
given by 



and 



bj{xi) = b T(0 + 1) L?{ Xi ), j = Q,...,K (110) 

aj { Xi ) = 0, j = 0,...,k-l (111) 

^•(^ = 7^^ + 1)^(4 j = k,...,n + k (112) 

with Xi being the solution of L^ +1 (x) — and 



4x 



^ 2 - a 2 ' 



(113) 
□ 



B.2 Theorem 1 



Now, we will use lemma [3] to get a solution of a set of m coupled recursion 
relations. 

Theorem 1. Consider a set ofm + 1 cut recursion relations for the generalized 
Laguerre polynomials with mixing terms, Xi> as described below, 



- (2j + a + 1 - x)a° j {x) + (j + l)(j + a + l)a$ +1 (a;) - Xo^+^x) = 0, 
J_ x (a;) - (2j + «i + 1 - a:)o}(a;) + (j + + ai + l)o) +1 (a?) - Xia' +9l (s) = 0, 
a^_ x (x) - (2j + a 2 + 1 - x)a 2 (x) + (j + l)(j + a 2 + l)a? +1 (ai) - x 2 af +g2 (a;) = 0, 



3 

a 



(x) ~ (2j + a m _! + 1 - x)aj-\x)- 

,m-l 



+ (.? + 1)0" + ttm-i + l)a™ /(.t) - Xm-iaf+fc^ (s) = 0, 
a™ ^z) - (2j + a m + 1 - ir)^) + (j + l)(j + a m + 1)0^(1) = 0, 
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Then, assuming that there are k + 1 coefficients aj l (x), there will exist K + l 
nontrivial solutions for x such that = 0, provided that for all integers 

< i < m we have, (k m = 0), 

Qm—i 

Proof. We start by solving the set of equations for a" l (x) using the results of 
lemma [T] Two situations can be possible, either a; is a solution to the equation 
= 0, or it is not. 

. L<% x {x)^Q 



— h 



(114) 



If L"'^ 1 (x) 7^ 0, lemma [T] states that all variables aj l (x) must vanish. In 
this case the equations for the variables a™ _1 (x) become homogeneous. 
Hence, we recover the assumptions of the present theorem with m replaced 
by m — 1. 



Lemma Q] predicts k + 1 nonzero solutions for all variables a™ (x) . 

In the following we assume without loss of generality the second possibility, 
x being a solution of L"^ 1 (x) = 0. The set of equations for the variables a™(x) 
and aj i ~ 1 (x) satisfies the assumptions of lemma [31 provided that 

1 (a m - a m _i) = /c m _i with k m _i a positive integer, 

(115) 



9m— 1 — 1 — km—1- 

According to lemma[3]we need k+2 — q m -i coefficients a™ -1 (a;) for consistency. 
Moreover, this theorem gives the solutions for a" l ^ 1 {x) as being proportional to 
aj l (x) and introduces the mixing constant 

4Xm— 1 /,-,»! 

7m- 1 = —o 2 ' ( 116 ) 



Q 



m — 1 



namely, 



a™-\x) = a™ 2 4 ^ m ' 2 1 r(a m + l)L^ fcm _». (117) 



Consequently, considering the set of equations for a™ 2 {x) we have 
afS x 2 (x) - (2j + a m _ 2 + 1 - x)a™- 2 (x) + 



3 

+ (j + + a m _ 2 + l)c%T*(x) + Xm-ao^_ a (a:) = 0. (118) 
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This can be solved using lemma [5] We assume that 



a™-\x) = 7m _ 2 a™ r(a m + l)a££_ 3 (a:) (119) 
and then, provided that, 

-(a m - a m _ 2 ) = k m -2 with k m -2 a positive integer, 

9m- 2 = 1 — k m -2 + k m -i. 

we find another mixing constant, namely, 

7m-2 = 2 _ 2 ■ (121) 

a m a m-2 

Again, we need k + 3 — <7 m _i — Qm-2 coefficients aj i ~ 2 (x). Hence, finally, 

a™- 2 (x) = a™ 7m _ l7m _ 2 r(a m + l)L^ fcm „», (122) 

which we write as 

a] n - 2 {x) = a^r m _ 1;m _ 2 r(a m + l)L°^ m _ 2 0r), (123) 

with 

a; 

r«,, = n7t. (124) 

t=y 

One can repeat these steps m — 2 times and find the solution for all a Ax). 

Summarizing, for all n + 1 roots of the equation L"™ 1 (cc) = 0, denoted by 
Xi, we have a nontrivial solution of the form, 

ay{xi) = a% T(a m + l)L^(x i ), 

with < j < k 
a™- 1 ^) = a™ r m _ 1>m _a I> m + l)Lp km i ( Xl ), 

with fc m _i < j < k + fe m _i 
af- 2 { Xl ) = a™ r ro _ 1>ro _ 2 r(a m + l)Lp km i _ km 2 (x t ), 

with fc m _i + fc m _ 2 < j < « + k m -i + k m ~ 2 



a](xi) = a™ r m _ X) i T(a m + l)££" E m-i fc( 

m—1 m— 1 

With ( ^ fc t ) < j < K + ( 51 

t=l t=l 
a% Xl ) = a™ r m _ li0 T(a m + 1)L Q "* ^ (a*). 

m — 1 m— 1 

With ( ^ fc t ) < j < K + ( X] k 



t=0 



□ 
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B.3 Application of theorem 1 to the recurrence in the 
bosonic sector 

We apply theorem Q] to the recurrence relation eq. (fT3|) . 
B.3.1 Even solutions 

In the case of even solutions we get the identification 



a t = 6t + 3, 

Xt = -f(t+l)(2t+l), 

h = 5 (at - a t _i) = 3, 

qt = -2 = 1- k t . 



Therefore, we obtain 



1 (t+l)(2t + l) 
7 * = "3(2 K +l)^-(2t + l)^ (126) 



a t = 6t + 6, 

X* = -f(t+l)(2t + 3), 

fct = |(a t - a t _i) = 3, 

Qt = -2 = l-fct, 



and obtain 



Thus, 



i (t + i)(t + j) 

7f 6 + l) 2 - (t + l) 2 ' 



(125) 



and 

B.3. 2 Odd solutions 

In the case of odd solutions, we can identify as follows 



(128) 



B.4 Theorem 2 

Eventually, we need the solutions to a twofold coupled Laguerre recursion rela- 
tions. We present it in a form of the following theorem. 



33 



Theorem 2. Consider a set ofm + 1 cut recursion relations for the generalized 
Laguerre polynomials with mixing terms as described below, 

S^-a](x) 
S? a ■ a){x) 



Xo 




- qo ( x ) 


- Mo 




= o, 


Xi 






- Ml 




= 0, 


X2 




-*» 


- M2 









i-3 




3 (x) 


Xm — 3 


a m ~ 2 


-,(») 


Pm — 3 "-j_|_p m 




= 




i-2 


<" 


2 (x) 


- Xm-2 




_,(*) 


., „m 
A*m— 2 u j+p m 


_,(*) 


= 




l — 1 






— Xm-1 


m 
a j+q m 


_» 


= o, 










- • a; 

















(131) 

// x is one of £/ie k + 1 roots o/ £/ie equation L"Ti (x) = i/ien, i/iere exzsi 
one nontrivial solution with variables a\(x) partially nonzero (for eachp at least 
for one value of the index t af(x) 7^ 0) , provided that for all integers < i < m 
we have, (kt = 0, for t > m) 

^(^m ^m—i) — ^m— ii 

_ , _ , , , (132) 



Qm—i 1 km—i ~t~ km- 



i+1- 



Proof. The prove goes in a similar manner as the proof of theoremQ] We start by 
solving the set of equations for aj 1 using the results of lemma [TJ By assumption 
we have that L"™ 1 (a;) = 0. The set of equations for aj l (x) and a 7 J l ~ 1 (x) satisfies 
the assumptions of lemma |3l provided that 

i(a m - a m _i) = k m - X with fe m _x a positive integer, 

9m— 1 — 1 — km— 1- 

Therefore, lemma [3] gives the solutions for aj l ' 1 (x) as being proportional to 
aj l (x) and introduces a mixing constant, 

4Xm— 1 noA\ 

7m- 1 = — 2 ' ( 134 ) 

m— 1 

Consequently, considering the set of equations for a™~ 2 (;r) we have, 

g a m _ 2 . a ,n-2 {x) _ Xm _ a0 ^_ a (x) - M m - 2 af +Pm _ 2 (s) = 0, (135) 

where 

o™- 1 ^) = a m 7 m-ir(a m + l)L^ fcm _», (136) 

o7*(i) = o% T(a m + l)L^(x). (137) 
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Thus, we have 

a;""' • a™- 2 (x) - X m-2lm-ia^T(a m + l)L^ km _ 1+9m _ 2 (x) 

- /x m _ 2 a£T(a m + 1)L"+ Pm _ 2 (a) = 0. (138) 

Provided that, 

~{otm — a ra -a) = A: m _ 2 with f; m _ 2 a positive integer, 

Pm-2 = 1 - fem-2, '- 139 "' 
<7m-2 = 1 — + »m— 1) 

we can use lemma [2] and obtain the solution for dj l ~ 2 (x) as being proportional 
to aj l (x) with a mixing constant given by, 

7m-2 = 5 5 ■ (140) 

cv — cv 

^ rn m— 2 

Next, the set of equations for a™ - (#) reads, 

S^" 3 ■ a-" 3 (x) - x«-3^t s W " Vrn-zaJ-^x) = 0, (141) 

where a™ -1 (2) and a 1 J l ~ 2 (x) are given now by 

a"- 2 (x) = a™ 7m _ 2 r(a m + l)^ m _ 2 (x), (142) 
af-^x) = a™ 7m _ x I> m + l)^^ (s). (143) 

Thus, we have 

go,.-, . a rn-3 {x) _ a ™ Xm _ 37m _ 2 r (am + l)Lp km _ 2+qm 3 (x) 

- a^Mm-37m-ir(a m + 1) J£" fcm _ 1+Pm _ 3 (*) = 0, (144) 

Again, provided that 



-(a m - a TO _ 3 ) = /c m _ 3 with fc m _ 3 a positive integer, 

Pm-3 = 1 — + 



(145) 



we can use lemma [5] and obtain the solution for a J 1 3 (x) as being proportional 
to a™ (a;) with a mixing constant given by, 

4(Xm-37m-2 + Mm-37m-l) ^.m 
7m-3 = 5 2 • ( 146 ) 

at — or q 

One can repeat these steps m — 3 times and find the solution for all a p -. 
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Summarizing, for all roots Xi of the equation L"^p 1 (x) = we have a non- 
trivial solution of the form, 

a? = T(a m + 1) L" m (xi), 

with < j < k 
af-^Xi) = a™ r m _ lim _! T(a m + ^Lp^ Jx,), 

with fc m _i < j < K + fc m _i 

af^ 2 (xi) = a™ r m _i, m _ 2 r(a m + i)L"™ kml _ km2 

with fe m _i + fe m _ 2 < j < k + 



a)(xi) = a™ r m _ M r(a m + l)L a ™ _ x (^), 



771— 1 m — 2 

z, 

t=l ' ' t=l 



with (J2 k t)<j< K +(J2 k 



oPAxi) = a% r m _ 1)0 r(a m + 1)L Q ™ (a*), 

3 2^t=o ft * 



771 — i m— ± 

with ( fc t ) < j : < K + ( ^ h) 



where we have introduced by force a similar notation as used in the previous 
theorem, namely, 

r*,, = 7», (147) 

with 7 y defined recursively by 

Xyly + l+»y7y + ^ ^ 
a — or 

and Jy> m = 1. □ 

B.5 Application of theorem 2 to the recurrence in the 
sector np = 1 

In this subsection we use the theorem [2] to find solutions of the recurrence 
relation in the sector with rip = 1. 

B.5.1 Even solutions 

In order to obtain solutions of eas.([6"5j). we identify as follows 

3f + 4, t odd, 
3t + 5, £ even, 

-|(t+l), Tiodd, (149) 



at = 



XL 



— 4(t + 1), t even, 



Mt = -|(t+l)(t + 2), 



36 



There is neither closed expression for j y nor for T x y . 
B.5.2 Odd solutions 

In a similar manner, solutions of eqs. (|64p . can be obtained after the following 
identification, 



at = 



3t + 4, t even, 
3f + 5, t odd, 



t 



_ t even, (150) 

Xt ~ \ -§(*+ 1), t odd, 
IH = -f(t+l)(t + 2). 

Similarly, in this case closed expressions for 7 y and for T x ,y do not exist. 
B.6 Corollary 

Corollary 1. The set of all possible values of the x parameter for which there 
exist a nontrivial solution of the cut set of generalized Laguerre recursion rela- 
tions does not depend on the precise form of the mixing coefficients x and n as 
long as \ and [i do not depend on j . 

Proof. The set of m + 1 cut generalized Laguerre recursion relations of the 
form presented, in theorems [JJ or [2J have a specific structure, namely, the i- 
th recursion relation contains mixing terms proportional only to the variables 
of the j-th recursion relation where j > i. Let us assume that Vj > i all 
the amplitudes described by the j-th recursion relation vanish. Then, the i- 
th recursion relation has no mixing terms and can be solved by the lemma [TJ 
Hence, it follows from the lemma [T] that the possible values of the x parameter 
are given by the equation 

X#+i(aO = 0, (151) 

with some appropriate ft,. Let us now assume that x is not a solution of the 
above equation. Then, the amplitudes described by the i-th recursion relation 
must all vanish. In this case, the (i — l)-th recursion relation has no mixing 
terms and can be solved by the lemma [TJ The set of possible values of the x 
parameter is therefore given by the sum of zeros of the equations, 

L%+i(x)=0, (152) 
L a K *Zl +1 (x)=0. (153) 

By induction, we can conclude that the set of all possible values of the x param- 
eter for which there exist a nontrivial solution of the cut set of m + 1 generalized 
Laguerre recursion relations is given by the equation 



(154) 
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We complete the proof by noting that the cut set of to + 1 generalized Laguerre 
recursion relations by assumption contains one recursion relation without mixing 
terms which can be solved with the use of lemma [T] □ 



